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A hydrodynamical description for vortex states in type II superconductors with no 
pinning is presented based on the time-dependent Gintzburg Landau equation (TDGL). In 
contrast to the familiar extension of a single vortex dynamics based on the force balance, 
our description is consistent with the known rotating superfluid hydrodynamics and, at 
long distance limit, with that following from the Landau level approach for high field case, 
and enables one to study nonlocal conductivities perpendicular to the field in terms of 
Kubo formula. Typically, the nonlocal conductivities deviate from the usual vortex flow 
expressions as the nonlocality parallel to the field becomes weaker than that perpendicular 
to the field measuring a degree of positional correlations, and, for instance, the dc Hall 
conductivity nonlocal only in directions perpendicular to the field becomes zero in the 
vortex lattice with infinite shear viscosity. Various situations are discussed based on the 
resulting expressions. 
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1. Introduction 

At present, understanding theoretically the nonlocality [1,2] of linear resistivity in a 
clean system (with no pinning) seems to be one of central problems on the vortex states of 
a type II superconductor. By analogy to the usual viscous fluid, Marchetti and Nelson [1] 
have previously proposed within a linear hydrodynamics that, deep in the liquid regime, 
a shear-viscous force proportional to — d\dt& will take the place of the shear term in the 
elastic force f e i in the following force-balance equation for the vortex lattice 

d t (-T 1 s + T 2 (z x s)) + J x B = f e i, (1) 

and that consequently, a positive nonlocal term proportional to q\ should exist in the 
conductivity perpendicular to B. Here s denotes the displacement field of vortex positions, 
J the external current, B the uniform flux density parallel to z axis, q± the wavevector 
perpendicular to B, and Ti and T 2 will be given later. This equation is essentially an 
extension [3] of the single vortex dynamics to interacting vortex states. Subsequently, 
the presence of other nonlocal terms was proposed [2] even in the case with no sample 
disorder on the basis of experiments in heavily twinned samples [4], and, deep in the liquid 
regime, the presence of a large conductivity nonlocal in the field direction was argued. It is 
important to theoretically clarify to what degree these proposals can be justified beyond the 
phenomenology [1,2,5] and within a basic dynamical model for type II superconductors such 
as the time dependent Gintzburg-Landau equation(TDGL) which reasonably describes the 
uniform linear dissipation in terms of Kubo formula [6]. Actually, the eq.(l) is not correct in 
the following senses. Firstly, it cannot be used in calculations of conductivity based on the 
use of Kubo formula [5,6]. Secondly, it is not clear from (1), in which f e i is independent 
of J and of the time derivative dt, how the possible viscous (nonlocal) terms in the dc 
conductivity are changed as the vortex lattice is formed. Further, the dispersion ~ of 
the shear mode found as a phase fluctuation formally and generally [7,8] cannot be seen 
in eq.(l). 

In the present paper, a correct hydrodynamical approach, or equivalently, a dynamical 
harmonic analysis on the transport in a vortex state, particularly the vortex lattice, is 
presented on the basis of TDGL. Our approach has no difficulties mentioned above and 
makes it possible to study nonlocal conductivities above and below the melting transition 
with no pinning, and several consequences on the resulting conductivities are discussed. 
As has been understood through studies [9] on nonlocal linear and static responses, the 
use of the harmonic analysis for 3D and layered systems often lead to misunderstandings 
on physics in the liquid regime, and hence we will not consider situations with negligibly 
small positional correlations of vortices. 

2. TDGL hydrodynamics 

Our starting point is TDGL describing dynamics of the order parameter ip 

a( 7 + W)d t ^ + m\ 2 - po) + ae (-id - —A) 2 ]i(j = 0, (2) 

<Po 

where curlA = Bz, 7, a, and b are positive, (j)Q is the flux quantum, £ the coherence length, 
and po is the spatial average of the mean squared order parameter. Since the linearized 
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version of (2) is considered throughout the paper, the noise term introducing fluctuations 
in equilibrium was neglected, and growths of time scales will be phenomenologically in- 
cluded later as the only fluctuation effects. Within the harmonic analysis given below, 
spatial variations of the gauge field have to be neglected in a consistent sense with the 
corresponding derivation of the infinite diamagnetic susceptibility in the Meissner state. 

Particularly upon calculations of linear responses, it is convenient to work rather 
within the harmonic approximation of the corresponding 'quantum' action with the gauge 
field disturbance 8 A (and with h = 1 ) 

S har = /[/3 5>( 7 |Q|+i 7 'Q)|^| 2 + [[ae \(-id-^(A + 5A(r)))ij(r)\ 2 

Jr V J r ^ 

+6(2|^ | 2 |^(r)| 2 + ^ *V(r) + c.c.))]]. (3) 

Since we are interested only in frequencies and wavevectors accompanying 5 A and sum- 
mations with respect them are not performed, examinations based on (3) are equivalent 
to the linearized analysis of (1) and formally independent of temperature. Here r is the 
imaginary time, i[)q, is the temporal Fourier transform of ip(r), Q is Matsubara frequency, 
j3 the inverse temperature, and the mean field solution of The superconducting (and 
real) part cr^ of dc linear conductivity tensor should be always calculated in terms of Kubo 
formula 



^j(q) 



<9p s (q, ifi ) 



<9fi 



d S^-p-HnTr^ exp(-S har )) 



n„^+o dfi SAi(q, Qo)SAj(-q, -fi ) 



(4) 

fi — +0 



where Qq is the external (Matsubara) frequency, i, j = x or y, and the absence of dc 
uniform (q = 0) superfluid density was used(see below). Since the ac and uniform case 
of (4) always becomes the usual vortex flow expression in the present formulation, the dc 
conductivity will be only considered hereafter. Since the harmonic modes appearing in (4) 
are spatial and temporal variations around the mean field solution linearly excited by the 
gauge field disturbance, the linear response quantities in vortex states resulting from such 
a harmonic analysis [10] are not accompanied by the temperature implying that this 
approach is not applicable to the linear dissipation parallel to £?, which is a consequence 
of critical fluctuations [6] . Therefore, considerations are largely focused on the vortex flow 
configuration. 

Before presenting general results, it is useful to first see results on the linear conductiv- 
ity perpendicular to B in the Landau level (LL) approach within the lowest and next-lowest 
LLs (see sec. 4 in Ref.7). In such a high field approximation the mean field solution and 
the shear elastic mode are found within the lowest LL, and the next-lowest LL, which 
does not participate in constructing the mean field solution [7, 9], provides the uniform 
displacement of vortices. Since, in constructing basis functions of eigenmodes, Matsubara 
frequencies in the 'quantum' action play similar roles to the wavevector parallel to B , 
including dynamical terms is easily performed following Ref.7. For later convenience an 
amplitude-dominated mode in the lowest LL, denoted by Sp, will also be included. Then, 
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to the lowest order in q±, the action (3) becomes S^l = Sq + Si, where 



S ~(3aY^ [ (TPoNlxd 2 + YM5p* Q xn + ^-\Spsi\ 
n Jr a 



+Poe o \d zX n - ^M 2 d 2 + ^r||aixd 2 ), (5) 

00 2a 

Si ^(3wj2 [ (Mn\\son\ 2 + l'n(s* m xs on ) z +e o \d z s Qn \ 2 

' ^'son + B-^Axonxz)! 2 ), (5)' 
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where s T (the transverse component of s) is given by r B (dx x -2) with vortex spacing r# = 
a/ 4>o/2tyB, Sq the uniform displacement with vanishing q±, \ a longitudinal phase variable, 
and 6*66 the resulting shear modulus. Variational equations with replacement ifi — > cj give 
eigenmodes of the vortex lattice. Note that the resulting dispersion a; (real frequency) = 
— i(a7 po) _1 C66 (Q-L^b) 4 for the shear mode[7,8] is different from that following from (1), 
and hence that the extrapolation [1,3] of the uniform (or equivalently, the single-) vortex 
dynamics to the shear mode is not justified within TDGL. Through the London limit (7) 
given below, the action (5) with the constraint s T = r 2 B {px x z) is found to correctly give 
hydrodynamical results at longer distances than 



I = r B \fc 66 r 2 B /2ap ^, 



which is typically of the order tb- 

The absence of the dc uniform superfluid density (helicity modulus) is obvious from 
(4) and (5)', implying that a constant twist pitch 6Aq is transmuted into an uniform 
displacement of vortices keeping the free energy invariant. That is, situation is different 
from an elastic matter with free energy invariant with respect to uniform displacements. 
From (4) and (5)', we easily obtain the mean field expressions on diagonal (i = j) and Hall 
(iy^j) vortex flow conductivities 

rr la n \ n ^ ^ &j + 7* g»j) / R N 

= 0, = PsO WQ{1 + rlq y 2r (6) 

where p s o = 2(27r£ /</ , o) 2 Po a is the mean field superfluid density in zero field. Due to the 
constraint mentioned above, (5) does not give any nonlocal (and longitudinal) corrections 
to (6) associated with the freezing to the vortex lattice. Consistently with (6), the nonlocal 
superfluid density, namely the real part of the linear response function p s (q, uj + i0 + ), 
becomes p s ( n ± = 0, q\\) ~ p 8 or B q 2 /2 in dc limit. As pointed out in Ref.9, this q 2 behavior 
does not change through a 3D melting transition [11]. 

Next let us here comment on eigenmodes of (3) in 2D and nondissipative (7 = 0) case, 
corresponding to the rotating superfluid 4 He at zero temperature described according to 
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the Gross-Pitaevskii equation. Consistently with Ref.12 and 13, the order parameter in 
(3) will be divided into the amplitude p and phase ip; ip = y/p exp(i</?) and, just as in 
the usual London approximation, the presence of the field-induced vortices will be taken 
into account through the topological condition[14] on <p = (dip, d T <p) by neglecting any 
fluctuation-induced vortices. As a result, we obtain the following harmonic action 

S P h = aj^ {ii8pd TX + p e (dx - r~ B \z x s)) 2 + ^-5p 2 + (^J) 2 ). (7) 

The only difference of this action from (5) is that the longitudinal current 

v s = 2ae (d X ~ r~ B \z x s T )) = 2a£ 2 v (8) 

is nonzero in (7), which is necessary in considering longitudinal linear responses nonlocal in 
directions perpendicular to the field. When q±l > 1, the minimal coupling in (8) between 
X and s T is removed, and in principle the phase fluctuation behaves like in zero field case 
at such short lengths. Variational equations resulting from (7) agree with those following 
from an 'em' analogy in Ref.12, where an importance of v s at nonzero frequency was 
stressed in a context of superfluid 4 He. The action (7) does not include the term 

Sinc = i / 7"S(sxH> (9) 

Jr,T Zr B 

which is necessary in recovering the limit of incompressible fluid (i.e., Sp = and v s = 0). 
This addtional term with the coefficient 7", which should become 7', cannot be detected 
in this phase-only analysis for (3) (Hereafter, 7" will be assumed to be equivalent to 
7'). The eigenvalues following from Sp^ + Si nc precisely coincide, when 7' = 1, with 
those derived by Sonin[13] and give dispersions u 2 h ~ (07') 2 bC§§{q± rg) 4 / (l + (q± /i nc ) 2 ) 
and u 2 ~ (2^q/7V 2 3 ) 2 (1 + ((/j_ /i nc ) 2 ) of two modes corresponding to the shear (massless) 
and compression (massive) elastic modes, respectively, where /i nc = ^bpo/2ar B /^o- The 
quadratic dispersion of o; s h at low qj_ is an origin of the destruction [7, 8] of the long- 
ranged phase coherence in the vortex lattice. The variational equation, with replacement 
—ir — > t(real time), with respect to the shear displacement of S'ph + S- 1IIC becomes the 
transverse part of (1) with no external current if Ti = 0, 1?2 = apQ^"2iiB / 4>q, and v s = 0. 
In this case the constraint v s = does not imply that in (5)' but rather corresponds to the 
limit of incompressible fluid. Following [13], this limit is valid only at shorter lengths than 
/inc, which is at most of the order of several vortex spacings in a type II superconductor 
near the melting line in a moderate field. 

The above agreement with the well-known rotating superfluid hydrodynamics justifies 
the presence of the minimal-coupling in v s between the phase field x and shear displacement 
s T . Based on this, we assume below that the v 2 term appearing in London limit(7) will 
be found at the static level by summing up many higher LLs in GL harmonic analysis 
around the mean field state. When combining this with (5) and (5)', we are naturally led 
to invoking the following action [15] appropriate to examining nonlocal conductivities: 

S* = / [Poa|0|( 7 \ Xq n\ 2 +7,^fW ) + Q 2 |s^| 2 ] 
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Jr PO 



X Z) 



+ r -4~\d z sn\ 2 + \d zX n - ^-SAatW (10) 
2 (po 

where 5A± (<5A^) denotes the longitudinal (transverse) part of the gauge disturbance de- 
fined within x-y plane. For later convenience the shear modulus is assumed to have possible 
frequency and wavevector dependences, and the time scale of s was phenomenologically 
changed taking account of a possibility that it may have nonlocal corrections ^ q — 7 > 0, 
due to an origin other than the freezing to the vortex lattice [1,2] and irrelevant to statics 
of vortex states at long distances. The amplitude mode leading to a 2 |x| 2 term was 
neglected. The term (9) has to be included when deriving a nonlocality of Hall conductiv- 
ity. Again, the variational equation with respect to s T of (10) is different from (1) with 
Ti = apo^nB /<f>o due to the presence of nonzero v s . Further, by neglecting the gauge 
disturbance 8 A and examining the dispersion of the shear mode in the vortex lattice, it is 
found that, due to the presence of the first term in (10), setting v s = from the outset, 
as in (1), always leads to an erroneous result on the dispersion. Rather, at low enough 
q, the constraint v s ~ is established, and consistently, the second term in (10) becomes 
unimportant compared to the first term. The first term in (10) having the same form 
as in the lowest LL case (5) is required within TDGL formalism and actually is justified 
because the x variable to be related to the shear displacement at small but nonzero q± 
was shown[7] to, up to the lowest order in q±, become a phase change around ifjQ even if 
higher LLs are fully included. 

3. Nonlocal Conductivities 

It is straightforward to, using (4) and (10), find the nonlocal conductivities. First, let 
us discuss the vortex lattice where Cqq (q = 0, O = 0) 7^ 0. In this case the longitudinal 
(II Q-l) P ar t of the diagonal conductivity and the Hall conductivity are given by 



r(s) r\ %4 + ^J 2 ^ + 2(lgJr B ) 

XX [Q) Ps °2e («|(l + « 2 r|/2) + iM) 5 



4 S » Pso^ (1 + ? J r | /2)(? 2 (1 + 9 2 r | /2) + ^ 2) - (11) 

using the length / oc y/C 66 defined previously. The transverse part (JxJ T of the diagonal 
conductivity merely becomes, even in the liquid regime, ^ q cr xx (q±_ = 0,(/||)/7 (see (6)), 
and hence, in contrast to (11), is not affected by the positional correlation [16]. We note 
that the term ~ q 2 + q\l 2 in denominators is the dispersion of the 3D shear (Goldstone) 
mode destroying the true off-diagonal long ranged order [7,8]. In general, the magnitudes 
of these conductivities significantly depend much on the relative size of \q\\ \ and \q±_\. For 

instance, when \q\\\ > l\q±\/rs ~ |<?j_|, both a x s J L and a x s J are well approximated by cr^ 
given in (6) with replacement 7 — > On the contrary, when \q\\\ < lqj_, a x s x L typically 
becomes the value 7p s o/<?i_£o a ^ zero Q\\ ( or m 2D case). This expression is the same as the 
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2D result in zero field [17] (see also [5]), reflecting the fact that it arises entirely from the 

phase variation x- Correspondingly, o~ x s J decreases like ~ (q\\/q±) 2 and vanishes at zero q\\ 
or in 2D case. In other words, as the nonlocality parallel to B becomes negligible compared 
to that perpendicular to £?, contributions of a finite 7 9 are covered and overcome by the 
perfect positional correlation, i.e., the presence of the Goldstone mode. This feature is 
lost, at fixed q± and q\\, as / becomes shorter, i.e., with increasing field, which is consistent 
with the absence of the nonlocal conductivity induced by the positional correlation within 
the lowest LL (see a sentence following (6)). In addition, note that, as seen in q» — > 
limit, the infinite shear viscosity of the vortex lattice does not mean a zero value of the 
nonlocal resistivity. We emphasize that the remarkable difference between the uniform 
(q± = 0) result (6) and this large conductivity in the case with zero q\\ and small q± is 
a consequence of the fact that, at zero q±, the shear displacement s T decouples with the 
zero mode x an d changes into the uniform displacement, and hence that it is theoretically 
correct. 

The above result that a xx ' ~ q ± and a xy ' = when q\\ = and q±_ is also 
applicable to a state with nonzero CQ${q ^ 0, |0| = 0) such as the 2D hexatic liquid state 
which may be possible in low enough fields (The 3D hexatic phase [1, 10] is inconsistent 
with the first order freezing to the vortex lattice and not possible in real superconductors, 
and hence need not be considered). It is valuable to examine these dc linear responses in 
real experiments in relation to determinations of the position of the freezing to the vortex 
lattice and of the existence or absence of the hexatic phase in 2D systems. 

When trying to understand the liquid regime (disordered state) in the present ap- 
proach, some comments are necessary. As pointed out in Ref.9, the harmonic analysis 
cannot be used in the situation where the contributions with nonzero reciprocal lattice 
vectors are negligible even if the amplitude fluctuation of the order parameter is negligible. 
In such situations the nonvanishing transverse diamagnetic susceptibility results from the 
static superconducting fluctuation and, in layered systems, shows a dimensional crossover 
due to a competition between the layer spacing and a finite phase coherence length which 
cannot be found in the harmonic analysis. Consistently, this dimensional crossover will 
be seen also in the nonlocal conductivity. On the other hand, it is practically difficult at 
present to understand possible nonlocal corrections, just above the freezing point to the 
vortex lattice, to the vortex flow expression (6) according to the fluctuation theory [6] . For 
these reasons, we will only consider the case just above the freezing point to the vortex 
lattice and with vanishing q\\ (see, however, the last section). 

In this case, the corresponding results to (11) are given by 

°xx (?±»?|| = 0) = p sQ - t 



2£o l + qlr%(% + m±r%/poa)/2>y' 



' S(?± ' = ° )= Ps °Wo l + alr%( % + m lr%/ Po a)/2Y (12) 

where the Maxwell form [18] CQs(q, ~ ??|0| with shear viscosity r\ for the dynamical 
shear modulus was assumed. Interestingly, the expressions (12) smoothly lead to (11) in 
zero gy for the vortex lattice by taking r\ to be infinity. The terms 1 + rj(q±r b) 4 Pqcl 
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in the denominators of (12) again originates from the spectrum ui ~ — iCQQ(q±rs) 4 of 
the shear mode, which cannot be found from (1). Consequently, at high q± of the order 
re -1 , both (11) and (12) are well approximated by (6), suggesting that, for such a rapid 
variation of the external current, not only collective effects but also the sample disorder 
is irrelevant because the time scale to be affected by the sample disorder (pinning) will 
be 7 9 and not be 7. In general, in the case with a pinning effect where ^ q will grow, the 
viscous effect accompanying the freezing transition becomes relatively negligible. Further, 
the expressions (12) suggest that the low temperature limits of the nonlocal conductivities 
for small q± in 2D systems with pinnings, where we have no transition at nonzero tem- 
peratures, again become the corresponding ( and above-mentioned ) results in the vortex 
lattice independent of 7 g . 

4. Comments and Conclusion 

Finally, we will discuss consequences of possible viscous effects [1, 2] in pinning- free 
systems, other than rj, which should appear in nonlocalities of 7 q . According to (11) 
and (12), when \q^\ <C \q±\, it is difficult to practically divide contributions of 7 g from 
the positional correlation, and hence, we will only consider the case with vanishing q± 
but nonzero q\\, where the diagonal conductivities can be always expressed by (6) with 
replacement 7 — > 7 g . Further we will focus on the 3D region below a dimensional crossover 
suggested in sec. 3 (We note that this dimensional crossover is quite different from that 
argued through experiments in Ref.4). For this case, it is often argued [19] that thermally 
activated cutting (and reconnection) processes among vortices induced by a nonuniform 
(Q\\ 0) current will lead to a very long length scale and significantly increase ~f q even if the 
(thermally-induced) entanglement is absent [2]. This picture suggests that the time scale 
grows unlimitedly even in the vortex lattice under cooling, and hence that the bent vortices 
cannot move at low enough temperatures. It should be noted that such an argument is also 
applicable to a nondissipative case (with zero 7 but nonzero 7') by imagining a nonlocal 
growth of 7'. However, it is even unclear to us if this is a correct argument in the context 
of the linear hydrodynamics. In superfluid 4 He[20] the reconnection process does not 
seem to need a remarkable slow dynamics. In addition, the system-size dependence of 
a temperature [4, 21] characteristic of apparently nonlocal vortex motions is inconsistent 
with the argument based on the thermally activated vortex cutting processes, because the 
measured size dependence inevitably means an algebraic decrease of the cutting barrier 
with increasing the system-size, although intuitively such a decrease of the barrier should 
not be expected. In relation to this issue, we emphasize that viscous terms appearing in 
7 9 cannot reduce to elastic (and static) terms in the vortex lattice and that the (if any) 
entanglement [1] making a growth of 7 9 possible must disappear deep in the vortex lattice 
[22]. In Ref.9, the presence of a large and positive nonlocal (~ q^) contribution in 7 g 
was questioned on the basis of the observation that, in contrast to the pinning-induced 
activation form of the time scale, any nonlocal (i.e., wavevector-dependent) growth of the 
time scale may be incompatible with the uniform vortex flow due to the nonlinearity in 
Kubo formula in the nonGaussian fluctuation theory [6]. Nevertheless, a possibility of a 
large 7 9 for nonzero q should be further searched in an extention of the fluctuation theory 
[6] which is not available at present. Experimentally, measurements in twin- free samples 
corresponding to those in [4, 21] should be performed. Actually, the data on the resistivity 
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parallel to B in [21] seem to intuitively contradict the data in a more 3D-like situation 
[23]. 

In conclusion, a hydrodynamics for vortex states in type II superconductors consis- 
tent with the rotating superfluid hydrodynamics has been presented in order to study the 
nonlocal conductivities perpendicular to the applied field. In particular, the vanishing Hall 
conductivity and divergent diagonal conductivity in the case with nonlocality in perpen- 
dicular directions to the field will be applicable to (if any) 2D hexatic state and useful in 
experimentally judging its existence or absence in type II superconductors. 
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